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ABSTRACT 


Context. The nature of dark matter (DM) is still under intense debate. Sub-galactic scales are particularly critical, as different, currently 
viable DM models make diverse predictions on the expected abundance and density profile of DM haloes on these scales. 

Aims. We investigate the ability of sub-galactic DM haloes to act as strong lenses on background compact sources, producing gravi- 
tational lensing events on milli-arcsecond scales (milli-lenses), for different DM models. For each DM scenario, we explore whether 
a sample of ~ 5000 distant sources is sufficient to detect at least one milli-lens. 

Methods. We develop a semi-analytical model to estimate the milli-lensing optical depth as a function of the source’s redshift for 
various DM models. We employ the Press-Schechter formalism, as well as results from recent N-body simulations to compute the 
halo mass function, taking into account the appropriate spherically averaged density profile of haloes for each DM model. We treat 
the lensing system as a point-mass lens and invoke the effective surface mass density threshold to calculate the fraction of a halo that 
acts as a gravitational lens. We study three classes of dark matter models: cold DM, warm DM, and self-interacting DM. 

Results. We find that haloes consisting of warm DM turn out to be optically thin for strong gravitational milli-lensing (zero ex- 
pected lensing events). CDM haloes may produce lensing events depending on the steepness of the concentration-mass relation. 
Self-interacting DM haloes can efficiently act as gravitational milli-lenses only if haloes experience gravothermal collapse, resulting 


in highly dense central cores. 


1. Introduction 


1 [astro-ph.CO] 27 Sep 2022 


> One of the most groundbreaking findings during the last century 
CX) was the discovery of a mass excess in nearby galaxies that could 
ON not be explained by the amount of ordinary (luminous) matter 
CY that was found to exist in those galactic systems. The observa- 
CY) tion of flat rotation curves in spiral galaxies (e.g., Rubin et al. 
— 1980; Bosma 1981; Corbelli & Salucci 2000), the discrepancy 
ON between observed velocity dispersion measurements and those 
© predicted by the virial theorem, in elliptical galaxies and globu- 
CN lar clusters (e. g., Zwicky 1933; Faber & Jackson 1976), and the 
N presence of collapsed structures at high redshift (e.g., Gunn & 
> Gott 1972) were smoking guns for the existence of a new ex- 
e= traordinary form of matter in the Universe, called dark matter 
(DM), which neither emits nor absorbs radiation. The DM hy- 
= pothesis has led to several correct predictions and has explained 
many observational discrepancies that had emerged in the past 
from the comparison of observational data with the first cosmo- 
logical scenarios that were dust-radiation-only models. Despite 
the great success of the DM model, its nature remains unknown, 
making it one of the most fundamental unsolved questions in 
physics. 


The most widely accepted scenario for the origin of DM is 
the so-called cold dark matter (CDM), a part of the standard 
ACDM cosmological model that has been remarkably success- 
ful in explaining the properties of a wide range of large-scale 
observations, including the accelerating expansion of the Uni- 
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verse (Perlmutter et al. 1999), the power spectrum of the Cosmic 
Microwave Background (CMB) (Page et al. 2003), and the ob- 
served abundances of different types of light nuclei (Cyburt et al. 
2016). However, the ACDM paradigm still presents some dis- 
crepancies with observations, mostly at small scales. Such small- 
scale challenges include, among others, the ’cusp/core’ problem, 
the missing satellites problem, the too-big-to-fail problem, and 
the angular momentum catastrophe (for a review see Bullock & 
Boylan-Kolchin 2017; see also Perivolaropoulos & Skara 2022). 


An appealing solution to those problems is to modify the in- 
trinsic properties of DM particles. During the past few years, nu- 
merous DM alternatives and ACDM extensions have been pro- 
posed by several authors, with the purpose to address some of 
the ACDM challenges. One of the most promising DM alterna- 
tives is the warm dark matter (WDM) model (e.g., Viel et al. 
2005; Lovell et al. 2012), where particles have rest mass on the 
order of a few keV, such as sterile neutrinos or thermal relics, 
that had non-negligible velocities at early times. Another very 
popular DM scenario is the self-interacting dark matter (SIDM) 
model where particles interact with each other (e.g., Spergel & 
Steinhardt 2000) having non-negligible cross-sections, of the or- 
der of ~ lcm?/g (e.g., Zavala et al. 2013). Other more exotic 
DM alternatives include; ultra-light axion dark matter (Schwabe 
et al. 2016), dark atoms (for a review see Cline 2021), and fuzzy 
dark matter (e.g., Kulkarni & Ostriker 2022). Although these lat- 
ter models are not examined here, the toolkit we have developed 
can be straight-forwardly adapted to any DM model for which 
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the redshift-dependent mass function and the density profile of 
haloes and sub-haloes can be calculated. 


The properties of the DM particle affect the formation of 
DM structures on all scales, their stability, as well as their evo- 
lution in time. In addition, the fundamental attributes of DM 
particles modify the primordial power spectrum describing the 
initial overdensity seeds of cosmological structures. So, differ- 
ences in the intrinsic DM particle properties between different 
models are expected to lead to measurable deviations in the re- 
sulting mass function of collapsed objects. For instance, models 
which include light particles, such as WDM, feature a sharp cut- 
off in the differential halo mass function below a critical mass 
scale, which depends on DM particle mass: the mildly relativis- 
tic velocities of WDM particles in the early Universe led to 
small-scale density fluctuations being washed out (free stream- 
ing) (e.g., Melott & Schramm 1985; Viel et al. 2005). The den- 
sity profile of DM haloes also turns out to be noticeably differ- 
ent from model to model. For example, virialized haloes made 
of WDM particles typically have lower central densities with re- 
spect to CDM haloes of the same mass, by virtue of their gener- 
ally later formation epochs (e.g., Lovell et al. 2012). 


The study of DM haloes below sub-galactic scales turns out 
to be particularly crucial for the exploration of the nature of DM. 
Nevertheless, it is extremely challenging to detect such haloes 
directly, in order to measure their number density in the Universe 
and/or examine their internal structure, since they might not even 
form galaxies due to their small size. So, the only possible way 
to explore and study them is through gravitational effects. 


One of the most promising methods of detecting sub-galactic 
DM haloes is strong gravitational lensing, where light that passes 
near a massive object (the lens) is being deflected, traveling a 
longer path than it would in the absence of the gravitational po- 
tential of the lens (e.g., Weinberg 1972). As aresult, when a com- 
pact background source (for instance, a radio loud quasar) emits 
radiation with the lens being in between the source and the ob- 
server and close enough to the line-of-sight, then the path of the 
light is affected strongly, resulting in the emergence of multiple 
images of the background source on sky with different magnifi- 
cations (e.g., Vegetti et al. 2012), provided its projected surface 
density exceeds a threshold. This effect is commonly known as 
strong lensing (see for example, Wright & Brainerd 2000). In the 
special case where the source displays intrinsic variability, ob- 
servable time delays between the different images (pulses) may 
occur (Zackrisson & Riehm 2010). 


Gravitational lensing can be used to detect compact objects 
(COs) that could not be detected otherwise, such as primordial 
black holes (PBHs) or dense DM haloes. Press & Gunn (1973) 
introduced the idea of assessing the cosmological abundance of 
COs through their strong gravitational lensing effect on distant 
background sources. They demonstrated that the cosmological 
mass density of COs can be constrained by deriving the fraction 
of lensed radio sources. Later on, Wilkinson et al. (2001) carried 
out a search for milli-lenses (gravitational-lensing images with 
milli-arcsec separations) in Very Long Baseline Interferometry 
(VLBI) observations of a sample of 300 compact radio sources, 
but no lensed systems in the mass range ~ 10°Mo to ~ 108 Mo 
were found. Their negative result allowed them to place an upper 
limit Oco $ 0.01 (95% confidence) on the cosmological density 
of COs in this mass range, concluding that the contribution of a 
primordial supermassive BHs population to the dark matter con- 
tent of the Universe is negligible. The currently on-going Search 
for MIlli-LEnses (SMILE) project (Casadio et al. 2021) expands 
the search for milli-lenses in the range ~ 10° Mo to ~ 10? Mo, 
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to a complete sample of ~ 5000 radio-loud sources using VLBI 
data. 

Motivated by the potential of the SMILE project, in this work 
we develop a novel method to exploit its upcoming results with 
the purpose to derive constraints on the nature of DM and dis- 
criminate between currently viable DM scenarios. Our approach 
is based on the concept of the lensing optical depth, represent- 
ing the probability for an observed source to be gravitationally 
lensed by a foreground mass distribution. The prescription for 
the implementation of this method can be found in Zackrisson & 
Riehm (2007). Recently, several authors have followed similar 
approaches to place limits on the abundance of primordial black 
holes, using Fast Radio Bursts (FRBs) (e.g., Leung et al. 2022; 
Zhou et al. 2022a,b), Gamma-ray Bursts (GRBs) (e.g., Kalantari 
et al. 2021), afterglows of GRBs (Gao et al. 2022) and compact 
radio sources (Zhou et al. 2022c). Here, we pursue the possibility 
of sub-galactic DM haloes acting as gravitational mill-lenses. 

We derive the expected number of milli-lenses in the source 
sample of the SMILE project for various DM models by calcu- 
lating the milli-lensing optical depth as a function of the source’s 
redshift. This in turn depends on the halo mass function, as well 
as on the projected surface mass density. Both of these physical 
quantities have noticeable differences between various scenarios, 
and hence the milli-lensing optical depth exhibits differences be- 
tween DM models. 

The layout of this paper is as follows. In Sect. 2 we describe 
our calculation of the milli-lensing optical depth. In Sect. 3 we 
discuss the analytic descriptions we use for the structure of DM 
haloes for various cosmological DM scenarios, and their corre- 
sponding mass functions. In Sect. 4 we present the results of our 
calculations, which we discuss in Sect. 5 


2. Lensing probabilities 


The principal result of any survey for lensing systems in the ob- 
servable Universe is the number of confirmed lensed images in 
a complete sample of sources (e.g., Myers et al. 2003; Browne 
et al. 2003). To maximize the constraining power of this product, 
we have to connect it to theoretical models that predict the expec- 
tation value of lensing events taking into account the differences 
in the abundance and density profile of DM haloes between vari- 
ous DM models. The most straightforward way to achieve this is 
to compute the lensing optical depth for any given DM scenario. 

The lensing optical depth depends strongly on the mass func- 
tion of gravitational lenses and on the surface density profile of 
each halo which essentially is related to the density profile. It 
also depends on the cosmology. In this paper, we fix the cosmo- 
logical parameters to be Hy = 100 h km s~! Mpc™!, h = 0.7, 
Qn = 0.3, Qa = 0.7, n = 0.97, 6-(0) = 1.674, and og = 0.8. The 
overall results, however, are not sensitive to small variations in 
these parameters. 


2.1. Milli-lensing optical depth 


In order for our results to be applicable to the SMILE project, 
we are interested in lenses that produce multiple images with an- 
gular separation on the order of milli-arcseconds (milli-lenses). 
Thus, we focus on lenses of masses (10°—10°) Mo. For the calcu- 
lation of the mill-lensing optical depth, we adopt the prescription 
of Zackrisson & Riehm (2007). We treat the lens as a massive 
object of mass M, with an angular Einstein radius 


GM, Dis 
(oa DoDos , 


Be = () 
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where Dos, Dıs, and Do; are the angular-diameter distances from 
the observer to the source, from the lens to the source, and from 
the observer to lens, respectively, with the lens being at redshift 
z while the source is located at redshift zs. Daz can be written as 


Dap(easz0) E D 
ka +X. =. ara) 

AB(ZA, ZB lra da HO 

where H(z) is the Hubble parameter, 

2 

H 

(=) = On(1 +2)? + On, 3) 
Ho 


with ©,,, Qa referring to the present values of the density pa- 
rameters for matter and dark energy, respectively. 

The milli-lensing optical depth for a source at redshift z, is 
given by 


T(Zs) = j 
0 


where o(M;, z, zs) is the lensing (effective) cross-section 


4nGM, DoD 15 
—— 5 
c? Dos ©) 


cdt 
dz 


10°M. 

dn(Mi, Z, Zs) 

dz f o(Mı, z, z) — dM, (4) 
106Mo ' dM i 


o(M), z, zs) = 16D? = 


and dn/dM_ is the differential lens mass function 


dn(Mı,z,zs) _ dM(M),z, Zs) dn(M, z) 
dM, 7 dM, dM ` 


In Eq. (6), dn/dM is the differential halo mass function (see 
§3.2). We caution the reader of the two different masses enter- 
ing Eq. (6): the lens mass M;, and the halo mass M. These two 
are not in general the same because only the part of the halo 
in which the projected surface mass density exceeds the critical 
strong lensing threshold can act as a gravitational lens. The crit- 
ical surface density value for a source at redshift zs undergoing 
strong gravitational lensing by a foreground DM halo (lens) at 
redshift z is 


(6) 


Mı = c? Dos 
o(Mı, Z, Zs) 4nG DoD is 


LerlZ, Zs) = z 10° Mo kpc™?. (7) 
Therefore, in order to calculate the halo mass M for given 
Mj, z, and zs, we demand a solution of the equation 


x(M; Mı, z) = Ler(Z, Zs), (8) 


where È is the projected halo surface density described exten- 
sively in Sect. 3. Solving this equation numerically, we obtain 
M(Mi,z, Zs). We use the central finite difference approximation 
to estimate the derivative dM/dM, which appears in Eq. (6). 

Given that the halo surface mass density is obtained after 
an integration of the density profile, it is clear that the results 
will differ significantly from model to model, since each DM 
scenario predicts a different density profile. 


2.2. Expectation value of lensing events 


Once we obtain the milli-lensing optical depth, we evaluate the 
expectation number of lensing events in the SMILE source sam- 
ple using 


Nsources 


$, 1- exp(-T@s;)). ©) 


i=] 


N, = 


For t < 1, we can approximate Eq. (9) by 


Ne) 
i=1 


T(Zs,i). (10) 


The source sample, as well as their corresponding redshifts are 
described next. 


2.3. SMILE sample 


The source sample considered in this study is the one of 
SMILE!: a complete sample built starting from the complete 
sample used in the Cosmic Lens All-Sky Survey (CLASS; My- 
ers et al. 2003; Browne et al. 2003), the most successful search 
to-date for gravitational lens systems at galactic scales using 
radio frequencies. The complete sample of 11685 sources pre- 
sented in CLASS is drawn from two other catalogs: the 5 GHz 
GB6 catalog (Gregory et al. 1996), and the 1.4 GHz NVSS cat- 
alog (Condon et al. 1998). The CLASS catalog contains sources 
from declination 0° to 75°, with a minimum flux density of 30 
mJy at 5 GHz, flat spectral index (< 0.5) between 1.4 and 5 GHz, 
and Galactic latitude (|b| > 10°). The complete sample of 11685 
sources has been initially followed up in CLASS with low res- 
olution Very Long Array (VLA) observations at 8 GHz. SMILE 
started from the complete sample in CLASS and selected sources 
with total flux density at 8 GHz > 50 mJy. The 4968 sources that 
satisfy such a requirement make a complete sample of flat spec- 
trum sources at declination [0°, +75°]. 

In order to obtain redshift measurements for sources in the 
SMILE sample, we used the Optical Characteristics of Astro- 
metric Radio Sources (OCARS) catalog (Malkin 2018), contain- 
ing redshift measurements of a large number of radio sources 
observed in different VLBI astrometry programs. Of the 4968 
sources in SMILE, 2781 have an optical counterpart within 3 
arcsec with redshift measurements in OCARS. For the remaining 
sources we searched for an optical counterpart within 3 arc sec- 
onds, with known redshift, in NED?. In total, we collected red- 
shifts for ~ 2/3 of sources in SMILE. For the remaining ~ 1/3, 
we generated redshift measurements randomly selecting values 
from the known redshift sample. Their distribution is shown in 
Fig. 1. 


3. DM haloes & Mass functions 
3.1. Halo size & structure 


The internal structure of dark matter haloes affects the lensing 
optical depth, since the threshold for strong gravitational lensing 
is associated with the projected surface mass density which in 
turn is related to the shape of the density profile. For a detailed 
review of the various mass densities see Zavala & Frenk (2019), 
while a thorough comparison among various density profiles can 
be found in Merritt et al. (2006) work. 

An important finding of the past decades is that spherically 
averaged DM density profiles in N-body cosmological simula- 
tions have a universal form (Navarro et al. 1997). Such density 
profiles are described by a simple functional form characterized 
by only two free parameters. The first one is the concentration 
parameter, denoted by ca, which quantifies how concentrated the 


! https://smilescience.info/ 

2 The NASA/IPAC Extragalactic Database (NED) is funded by the Na- 
tional Aeronautics and Space Administration and operated by the Cali- 
fornia Institute of Technology. 


Article number, page 3 of 11 


A&A proofs: manuscript no. main 


C] known redshifts 


DZ randomly generated 


Probability density 


Zs 


Fig. 1. Redshift distribution for sources used in this study. Cyan solid 
line represents the redshift distribution of sources with known redshift, 
whereas golden dashed line stands for the distribution of the randomly 
selected redshift measurements from the known redshift sample. 


mass is towards the center of the halo. The other one is the char- 

acteristic radius, r;, which determines the distance from the cen- 

ter above which the density profile becomes steeper, i.e. quanti- 

fies roughly the size of the core. These two parameters are related 

to each other through 
TA 


CA= > 
rs 


(11) 


where ra is the virial radius. 

Since the distribution of mass is continuous, the boundary of 
a halo cannot be defined precisely. So, another major challenge 
is to come up with a robust method of determining the size of 
a halo uniquely. Thus far, numerous papers that deal with this 
problem have been published by several authors (Cole & Lacey 
1996; White 2001; Cuesta et al. 2008; Zavala & Frenk 2019). 
In general, the radius of a halo can be defined through the over- 
density parameter, A(z), which in principle depends on the cos- 
mology (Bryan & Norman 1998; Tinker et al. 2008; Naderi et al. 
2015; Seppi et al. 2021). In particular, it represents the radius 
where the mean interior density is A(z) times the critical density 
of the Universe p,,(z), namely 


3 a 
i | Peer = Apo, (12) 
dar, 0 
where the critical density is given by 
3H(z) HOY 
cr = = Pero | y] > 13 
Perk?) = -og TP Ho (13) 


with G being the Newtonian gravitational constant, while pero 
accounts for the critical density of the Universe at redshift z = 0, 
and H(z) is given in Eq. (3). 

The halo mass M,, which is the mass contained within a 
sphere of radius ra, is given by 


4n 3 
Ma = A>PryPers 


14 
3 (14) 
and as a result the halo radius can also be written as 
3M, 1/3 
M, z) = | —— i 15 
ra(Ma,z) (y) (15) 
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However, the most commonly used way to determine the halo’s 
size is to consider that the overdensity parameter J is fixed and 
equal to 200, since it turns out to be a rather convenient way to 
define the boundary of the halo and simplifies the calculations 
(e.g., Cole & Lacey 1996). Taking this fact into account, we fix 
the overdensity to be 4 = 200, throughout this paper, and there- 
fore the halo mass is Myo (hereafter, M), the halo radius is r200, 
and the concentration is C299 (hereafter, c). 

Other quantities used extensively below are the enclosed 
mass, Menc, the projected surface mass density, Ł, and the lens 
mass, Mı. The enclosed mass is defined as the mass which is 
contained within a sphere of radius r, namely 


Mont) = 4 | Adr, (16) 
0 


where we assume that we deal with spherically symmetric ob- 
jects. The projected surface density is derived simply by the in- 
tegration of the mass density along the line of sight (Wright & 
Brainerd 2000; Mo et al. 2010; Dhar & Williams 2010; Retana- 
Montenegro et al. 2012; Lapi et al. 2012) 


(8) = f oaz, (17) 


where s = Vr? +z? is the projected radius (orthogonal to the 
line-of-sight) relative to the center. 

Having defined the projected surface density we can infer the 
gravitational lens mass M;, by carrying out an integration of the 
ÈŁ(s) over the disk which has radius s (see for example Eq. 41 in 
Retana-Montenegro et al. 2012) 


Ms) = 27 { EO dx. (18) 
0 


This quantity is the mass contained within an infinite cylinder 
of radius s in which the mass distribution is characterized by 
the mass density profile p(r). Eq. (8) cannot be solved indepen- 
dently, but must simultaneously satisfy Eq. (18), owing to the 
fact that one of the input parameters in Eq. (8) is the lens mass. 
Therefore, the relation between lens mass and halo mass that is 
needed in the computation of the lensing optical depth is derived 
only after solving this coupled non-linear system. 


3.2. Mass function 


The computation of the lensing optical depth depends on the 
number density of lensing objects which in our case are DM 
haloes. So, we need to obtain a formula that determines the dis- 
tribution of virialized dark matter haloes (in the field) per vol- 
ume element and per halo mass at a given redshift z, in order 
to estimate the lensing optical depth. This problem has been ad- 
dressed by several authors in the past, either analytically (Press 
& Schechter 1974; Bond et al. 1991; Pavlidou & Fields 2005) or 
numerically (Jenkins et al. 2001; Tinker et al. 2008). All these 
works are based on the spherical collapse scenario (e.g., Gunn & 
Gott 1972; Naderi et al. 2015). Improvements using ellipsoidal 
collapse do exist (e.g., Sheth et al. 2001), but we will not use 
them in this paper. 


3.2.1. CDM mass function 


The differential halo mass function of CDM haloes (the number 
of haloes with mass between the range M and M+dM per proper 
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volume at a given redshift) reads (Press & Schechter 1974, see 
also Appendix A) 


dN _ 
dM 


2 Pm) 6¢(Z) 


n(M, z) = 
(M, 2) nM A 


2 

aad ex - o| : (19) 
dM 20%4 

where M refers to the halo mass, z is the redshift, p,,(z) is the 
mean matter density of the Universe at redshift z, and 6,(z) de- 
notes the overdensity of a structure collapsing at redshift z lin- 
early extrapolated to the present. Moreover, oy is the rms of the 
density field smoothed on scale M (see Appendix A). 


3.2.2. WDM mass function 


Although the halo mass function for WDM haloes shows a simi- 
lar behavior to the CDM one on galaxy clusters scales, it exhibits 
a cutoff below the dwarf galaxy scale, owing to the free stream- 
ing of WDM particles in the early Universe. The most commonly 
applied method to derive the WDM halo mass function is the de- 
velopment of N-body simulations that evolve the primordial den- 
sity field perturbations in time, leading to collapsed DM haloes 
(e.g., Schneider et al. 2012; Bose et al. 2016; Lovell 2020a,b). In 
this study, we choose to use the numerical fit offered by Lovell 
(2020b) to take into account the cutoff in the mass function of 
WDM haloes with respect to the CDM one. The halo mass func- 
tion in the case of WDM is given by 


Mim B id 
nwou(M,2) = ncpu(M, ali + (n) ) , (20) 


where Mj, is a characteristic mass scale (the half-mode mass), 
while a, 8, y are parameters of the fit that have been found to 
be 2.3, 0.8, —1, respectively. The half-mode mass is associated 
with the free-streaming length which in turn is related to the rest 
mass of the WDM particle. 

Here, we are interested in exploring the case where the 
WDM is made of collision-less particles (thermal relics) of mass 
Mwpm = 3.3 keV. In this scenario, the theoretical value for the 
half-mode mass is Mam = 2 x 108 Mo (e.g., Bose et al. 2016). 
This value for the half-mode mass coincides with the one for the 
well-motivated sterile-neutrino model in which particles are as- 
sumed to have rest mass equal to 7 keV and lepton asymmetry 
number Le = 8.66 (see Bose et al. 2016). Sterile neutrinos are 
part of the neutrino Minimal Standard Model (vMSM; Boyarsky 
et al. 2009) which is a simple extension to the Standard Model 
of particle physics. It has been introduced to explain the uniden- 
tified 3.53 keV X-ray line observed recently in galaxies (see for 
example, Boyarsky et al. 2014) by considering this line to be the 
decay signal of those 7 keV sterile neutrinos. 

Given that the cutoff in the halo mass function of these two 
WDM models is determined by the same half-mode mass and 
that the internal structure of haloes is identical, the inferences 
of this work concerning the 3.3 keV thermal relic WDM par- 
ticle will also be valid for the 7 keV sterile neutrinos model. 
In Fig. 2, we display the differential halo mass function for the 
CDM model and for the WDM one investigated here as a func- 
tion of the halo mass for various redshifts. Solid lines correspond 
to the proper density of haloes of mass M at different epochs 
(redshifts) divided by the present critical density of the Uni- 
verse, whereas dashed lines represent the same quantities, but 
for WDM. A major difference between the CDM mass function 
and the mass function of WDM is that the latter one exhibits a 
cutoff at the dwarf-galaxy scales (~ 10° Mo). This distinctive 
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Fig. 2. Comparison of the CDM differential halo mass function with the 
one of the WDM model for various redshifts. Solid lines refer to CDM, 
while dashed lines correspond to WDM. The vertical axis corresponds 
to the halo mass while the y-axis shows the proper density of haloes of 
mass M normalized to the present value of the critical density of the 
Universe. 


feature of the halo mass function in the WDM scenario has a 
considerable implication to the ability of DM haloes to act effi- 
ciently as gravitational milli-lenses on background sources, be- 
cause the mass function is directly involved in the calculation of 
the milli-lensing optical depth (Eq. 4). 


3.3. CDM halo density profile 


In this study, we employ the Navarro-Frenk-White (NFW) pro- 
file for the description of the mass distribution within CDM 
haloes (Navarro et al. 1995, 1996, 1997) 


be 
Ir Err? 


where r, is the characteristic radius, while ¢, is calculated by 


P(T, Fs, C2) = Per(Z) (21) 


200 o 
3 In(l+c)-c/(1 +0)’ 


bc = (22) 
with c being the concentration parameter, which is not mass in- 
dependent, but correlates strongly with the halo mass, as well 
as with the redshift, following a simple scaling law (see for ex- 
ample, Bullock et al. 2001; Neto et al. 2007; Prada et al. 2012; 
Dutton & Maccid 2014; Klypin et al. 2016; Shan et al. 2017; 
Ragagnin et al. 2019, 2021). The NFW profile predicts a cuspy 
halo’s center, since the mass density goes as ~ r~! near the center 
of the halo. 

Using Eq. (21) along with Eq. (16), we obtain the enclosed 
mass 


M 
Mence (M, z, M,z)) = ————— Fir, c, , (23 
(M, z, r200(M, z)) nite cite (r,c, 1200), (23) 
where 
AEE RLL R 24) 
1+ cr [r200 


From Eqs. (17) and (21), we obtain the surface density 


PcPer (z) T200C— 


1 
ÈZ (s, M, z, r2o0(M, z)) = 25 z S(s, c, r2x(M,z)), (25) 
c*(s/r290)* — 1 
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where we have defined for convenience 


-1 
cos” (ra00/CS . 
pa OS wie) 200/¢5) if s >r, 


yc?(s/r200)} — 1 


SeroMD = V y (26) 
- _——._ ifs<r, 
vl —c?(s/r200)* 
From Eq. (18), the lens mass then is 
Mins (5, M, z,r200(M,2)) = H(s,c, r20) 
ens (S, M, Z, F > = Spas a S, C, F > 
l Z, r200(M, Z ESEE 200 
(27) 
where 
H(s, c, r290) = In(cs/2r200) 
2 k -1 
—— arctan esiTao9 = 1 if cs > r200 
V(cs/r290)? — 1 cs/r200 + 1 
1 ifcs = F200 
2 1- 
———- arctanh 1~ es/r200 if cs < r200 
V1 = (es/r200) es/T200 + 1 
(28) 


3.4. WDM halo density profile 


WDM is made of particles that had non-negligible thermal ve- 
locities at early times. This major difference is expected to have 
an impact on the concentration of mass near the center, but not 
in the shape of the distribution of mass within a halo. Indeed, the 
density profile in WDM models can be well described by a NFW 
profile (see, e.g., Lovell et al. 2014; Bose et al. 2016). However, 
DM haloes consisting of WDM are typically formed at smaller 
redshifts with respect to the formation of CDM haloes. This dif- 
ference affects the concentration of the halo, which generally re- 
flects the density of the Universe at the epoch of halo formation 
(for a detailed discussion see, Schneider et al. 2012; Bose et al. 
2016; Zavala & Frenk 2019). Therefore, we assume that the mass 
distribution in WDM haloes is consistent with the NFW profile, 
but more fuzzy, i.e., less concentrated around the center. In order 
to calculate the concentration parameter for WDM haloes as a 
function of the mass and the redshift we use the findings of Bose 
et al. (2016). They offer the following simple functional form for 
the concentration parameter 


M, “Ea 5 
CWDM _ (1 bai un) (1 + 2, 
CCDM M 


(29) 


where yı = 60, y2 = 0.17, and B(z) = 0.026z — 0.04. Mpm is the 
half-mode mass and in this work is set to be Mam = 2 x 108 Mo 
which corresponds either to the model of thermal relics WDM 
particles of rest mass mwpm = 3.3 keV or to the 7 keV sterile 
neutrinos model (an extension to the Standard model). Due to 
their smaller concentrations, WDM haloes will be less likely to 
exceed the lensing surface-density threshold, resulting in a lower 
milli-lensing optical depth. 


3.5. SIDM halo density profile 


SIDM was originally introduced by Spergel & Steinhardt (2000) 
to explain observations of central densities in galaxies within the 
Local Group. Since then, numerous authors have argued that the 
self-interaction of particles leads to a core-like profile rather than 
a cusp-like profile. Such a feature could alleviate the cusp/core 
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problem arising for CDM, and hence is considered a well moti- 
vated DM alternative. 

There are two kinds of SIDM theories. In the first case the 
scattering rate per particle, I(r), is velocity independent which 
implies that the ratio of the effective cross section, ø, to the dark 
matter particle’s mass, m, is constant (e.g., Rocha et al. 2013; El- 
bert et al. 2015). In the second scenario, the scattering rate is ve- 
locity dependent and falls rapidly as the velocity increases (e.g., 
Zavala et al. 2013). For a recent discussion on SIDM models, 
as well as on the observational constraints on the self-scattering 
cross section, see Tulin & Yu (2018). 

In this work, we will assume that the scattering rate per par- 
ticle is velocity-independent and has the following form 


T(r) œ p(r(o/m)vyms (1), 


where p(r) is the DM mass density at radius r, while v,,,5 is the 
rms speed of dark matter particles. We consider a typical value 
for the ratio o/m ~ 1 cm*/ g, since SIDM models with smaller 
values, on the order of 0.1 cm?/g, are very similar to the CDM 
models even on scales smaller than dwarf galaxies and cannot 
produce detectable deviations from CDM predictions (Zavala 
et al. 2013). On the other hand, higher values of the cross section 
per mass, ~ 10 cm?/g, have already been ruled out by cluster ob- 
servations (see, e.g., Dawson et al. 2012). 

For the structure of SIDM haloes, we again assume spherical 
symmetry, but now we use a core-like profile. In fact, the mass 
density is well approximated by the Burkert profile (see, Burkert 
1995) which has also two free parameters and is given by the 
following formula 


(30) 


Pb 
(1 +r/rp) (1 + (r/r) = 


where rp is the scale (core) radius, while p, is the central den- 
sity. As in the NFW profile, the free parameters of the Burk- 
ert profile scale with the halo mass. In the special case where 
a/m ~ 1 cm?/ g, Rocha et al. (2013) have provided a couple of 
simple scaling laws that connect both the r, and p, with the halo 
mass, using data from N-body simulations. These relations (Eq. 
17 and Eq. 20 in Rocha et al. 2013) are given below 


PB(T, Tb, Pb) = 


" mM.. \8 
4 2.21 a 2 
1 kpc (oc) , we 
-0.19 
Moir 
B o k ) (33) 
Mo/pe 10'° Mo 


Since these relations have been derived using the virial mass, 
Mir, instead of M299 which we have employed throughout this 
work, we have to re-scale the density profile to be consistent with 
Eq. (12). Eq. (31) can be recast as 


Pb 
, 34 
(1+ r/ry) (1 + (/7p)?) os 


r200\ 1l 500 -1{ 7200 
—]+ = Inf + =] - tan | — 
rp 2 r rp 


(35) 


pB(r, M) = A(M) 


where 


-1 


Mr,? 
A(M) = 5 in + 
TPb 


and now we can use the halo mass M = Moo instead of the 
virial mass in Eqs. (32), (33). The term A(M) has been derived 
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Fig. 3. Comparison of mass density profiles for a given halo of mass 
M = 108 Ms at redshift z = 0 for various dark matter scenarios. The 
concentration for the CDM case is given by Eq. (38). 


by requiring the mean density inside a sphere of radius r20 to be 
200p-,; (see Eq. 12). 
Using Eqs. (16) and (34), we obtain for the enclosed mass 


2 
3 r r =i r 
Menc(r) = A(M appr, hafi + 5) + 2m(ı + z) — 2 tan (z) 
r 


Fp Yb 


b 
(36) 


The surface density cannot in general be derived analytically and 
therefore we have to perform the integration numerically. In the 
special case where s = 0, i.e., for the column density through 
the line-of-sight, the integration that returns the surface density 
(0) yields the closed-form expression 
m 
Zz(s = 0, M) = A(M)> porn, (37) 
where the index B indicates that this surface density arises from 
the Burkert profile. The surface mass density is maximized when 
s = 0 since p(r) is a monotonically decreasing function of r, and 
as a result for a SIDM halo of given mass, the maximum value 
of the surface density is determined by Eq. (37). This feature 
is of great importance in strong gravitational lensing where the 
surface density must exceed a critical threshold in order to sig- 
nificantly bend a light ray. 

In order to have a qualitative picture of the differences be- 
tween the three mass density profiles mentioned above, in Fig. 
3, we apply them to a DM halo of mass M = 108 Mo at redshift 
z = 0. For SIDM particles the density profile near the center is 
flat, while for CDM particles the profile near the center goes as 
r`, since we have used the NFW profile. For WDM particles the 
profile is NFW-like but with larger characteristic radius than in 
the CDM reflecting the fact that the concentration in the WOM 
scenario is smaller than the one in CDM. 


3.5.1. SIDM core collapse 


Even though most of SIDM models are in favor of a less dense 
core-like halo center, the strong self-interaction developed be- 
tween particles in the innermost region of the halo might have 
important implications in the dynamical evolution of the halo. 
In one scenario, strong self-interactions between particles in- 
duce a negative heat capacity, eventually leading to the forma- 
tion of a dense central core in the inner part of the halo (see e.g., 


Yang & Yu 2021, 2022). Yang & Yu (2021) demonstrated that 
such a scenario can be successful in explaining the observational 
excess of small-scale gravitational lenses in galaxy clusters re- 
ported in Meneghetti et al. (2020). They exploited the fact that 
at late stages of the gravothermal evolution of a halo composed 
of SIDM, the core might undergo gravothermal collapse, result- 
ing in a highly dense halo center, thereby increasing its lensing 
effect on background sources compared to CDM haloes. 

In the most extreme case, the collapsed core can further con- 
tract, eventually leading to the formation of a supermassive black 
hole (SMBH) at the halo center. This scenario was firstly pro- 
posed and studied extensively by Feng et al. (2021) (see also 
Feng et al. 2022) as a possible mechanism to explain the exis- 
tence and origin of SMBHs at high redshifts (z ~ 6 — 7). SIDM 
offers a natural mechanism for triggering dynamical instability, 
a necessary condition to form a black hole. This scenario can be 
tested and well-constrained through milli-lensing, since the cen- 
tral SMBH can effectively act as a strong gravitational lens and 
produce multiple images of a compact background source. 

Given that studies dealing with the core collapse scenario 
do not provide an exact formula for the final mass distribution 
of DM inside the collapsed halo, we shall restrict ourselves in 
investigating here only the latter, most extreme, scenario of core 
collapse where the formation of a SMBH takes place from the 
gravothermal collapse of the core. The exploration of this model 
yields an upper limit on the expectation value of lensing events 
in the SMILE source sample in the case of the SIDM scenario. 


4. Results 
4.1. CDM 
4.1.1. CDM: Model A 


We start by investigating the CDM scenario using a 
concentration-mass relation derived from N-body simulations. 
We employ the relation given in Ragagnin et al. (2019) to de- 
termine the dependence of the concentration parameter c on the 
redshift z, as well as on the halo mass M: 


si = M -0.12 1.47 0.16 
eee SN 1 Mo +z) ` 


In Fig. 5 we plot with a blue solid line the milli-lensing opti- 
cal depth obtained for this c(M,z). The value of the milli-lensing 
optical depth is well below ~ 1074 implying that even a sam- 
ple of ten thousands distant (z ~ 5) compact sources is highly 
unlikely to produce at least one lensing event. Indeed, perform- 
ing the summation in Eq. (10) over all sources involved in the 
SMILE sample, we end up with the value (Nexp) = 1.5 x 107, 


which makes detection of a milli-lens improbable. 


(38) 


4.1.2. CDM: Model B 


As a limiting case of the possible effect of the concentration- 
mass relation on our results, we also test a power-law extrapo- 
lation to lower masses of the empirical (fitted from observations 
rather than simulations) c-M relation shown in Fig. 13 of Prada 
et al. (2012) 


1 
log c(M, z) = 4.23 — 0.25 log(M/Mo) — 0.16 log (H) . (39) 


Regarding the dependence on redshift, we consider that it is 
identical to Eq. (38), but we stress that most studies suggest a 
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Fig. 4. Concentration-mass relation at z = 0 for different dark matter 
scenarios. 


weak dependence of the concentration on redshift, so even if we 
slightly modify the last term in Eq. (39) associated with the red- 
shift dependence, the overall results do not change noticeably. 
In practice, the concentration parameter is set by the halo mass. 
We note that Eq. (39) predicts higher values of the concentra- 
tion parameter with respect to the ones inferred from N-body 
simulations. Although this c — M relation has been derived from 
galaxy cluster observations and might overestimate the c param- 
eter of haloes on sub-galactic scales, recently Sengiil & Dvorkin 
(2022) investigated the strong lens system JVAS B1938+666, 
concluding that sub-galactic DM haloes can be highly concen- 
trated (c ~ 60), in line with Eq. (39). 

Using this relation in Eq. (4), we obtain the green dash- 
dotted line in Fig. 5, showing the milli-lensing optical depth as 
a function of the source redshift. Subsequently, using Eq. (10) 
to compute the expectation value of lensing events in the source 
sample of SMILE, we obtain (Nexp) ~ 1.2. This value deviates 
remarkably from the one corresponding to model A (see Sect. 
4.1.1), demonstrating that the concentration-mass relation plays 
a crucial role in the process of strong gravitational milli-lensing 
and can thus be strongly constrained with milli-lensing observa- 
tions. This value also places an upper limit on the expectation 
number of detected milli-lenses in the SMILE’s source sample, 
in the case where the properties of DM particles are in line with 
the framework of the CDM model. 

A comparison between the two concentration-mass relations 
related to the CDM scenario for redshift z = 0 can be found 
in Fig. 4. The concentration-mass relation given by Eq. (38) 
(Model A) is displayed with a blue solid line, while the green 
dash-dotted line stands for the c(M) considered in model B (i.e., 
Eq. 39). 


4.2. SIDM 


In the SIDM model, there are two possibilities which lead to 
quite different internal structure of haloes. 

The standard scenario is the one where the inner part of 
haloes is characterized by a core-like profile yielding the pro- 
jected surface mass density of Eq. (37). Since the surface mass 
density is maximized at the center of the halo (as long as the 
density profile is a decreasing function of r), if the central re- 
gion does not exceed the critical threshold for strong lensing, 
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then the halo will not act as a strong lens. Using Eq. (37) and 
Eq. (7), we conclude that the halo mass of a SIDM halo must 
be > 10/4 Mo, for the surface mass density at the center to ex- 
ceed the strong lensing threshold. However, this mass scale cor- 
responds to galaxy clusters and therefore no trustful inferences 
can be done without including the effect of strong lensing due 
to the presence of baryons. The main finding is that SIDM-only 
sub-galactic haloes cannot produce milli-lensing images since 
they are not dense enough to satisfy the strong lensing criterion. 


4.2.1. SIDM core collapse: Model C 


The second scenario related to SIDM haloes is based on the 
gravothermal core collapse process that might take place in the 
inner parts of SIDM haloes (see Sect. 3.5.1). 

Assuming that the halo in the beginning was described by 
a NFW profile with the concentration-mass relation to given by 
Eq. (38), we can calculate the mass enclosed inside a projected 
disc with radius equal to the scale radius r,;. Then, we consider 
the most extreme case where the entire core collapses into a very 
small but extremely dense core that eventually results in the for- 
mation of a compact object. This collapsed core is the part of the 
halo that can produce strong gravitational lensing of light emit- 
ted by background sources. 

In Fig. 5, we show with the black dotted line the milli-lensing 
optical depth in the case of SIDM core collapse. Having obtained 
the optical depth, we carry out the sum shown in Eq. (10) over 
the redshifts of the SMILE project sources and find the value 


(Nexp) = 13. 


4.3. WDM: Model D 


In order to derive the milli-lensing optical depth for the sce- 
nario of WDM, where particles are supposed to have rest mass 
Mwpm = 3.3 keV (thermal relic) or be sterile neutrinos with rest 
mass equal to 7keV, we take into account that the halo mass 
function is different from that in CDM, and so we adopt the fit 
offered by Lovell (2020b) (Eq. 20). Regarding the density pro- 
file we again employ the NFW one, but with the concentration 
parameter to be given by Eq. (29). This concentration-mass re- 
lation is however a fit that relates the concentration of WDM 
haloes with the one corresponding to CDM haloes, so we simply 
consider that the concentration-mass relation of CDM haloes is 
the one shown in Eq. (38) and in such a way we obtain a formula 
for the concentration of WDM haloes as a function of the halo 
mass and redshift. In Fig. 4, we display the concentration-mass 
relation which corresponds to WDM haloes at z = 0 with a red 
dashed line. 

Performing the integration of Eq. (4), we find the milli- 
lensing optical depth in the case of WDM, shown in Fig. 5 with 
the red dashed line. Combining this result with Eq. (10), we com- 
pute the expectation number of detected WDM milli-lenses, ob- 


taining (Nexp) = 1.1 x 107°. It is therefore extremely unlikely 


to detect any milli-lenses with SMILE if DM is in the form of 
WDM. 


5. Discussion & Conclusions 


In this work we have explored the ability of sub-galactic DM 
haloes to act as milli-lenses on background sources resulting 
in multiple images of the same source with angular separation 
on the order of milli-arcseconds, considering different DM mod- 
els. We have developed a semi-analytical method to estimate the 
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Fig. 5. Lensing optical depth as a function of the source redshift for 
different dark matter scenarios. 


expectation value of detected milli-lenses in several DM sce- 
narios, computing the lensing optical depth. We have modeled 
the number density and internal structure of haloes using either 
(semi)analytical calculations or fits to N-body simulation results, 
depending on the DM model. We have restricted ourselves in ap- 
plying the point-mass lens approximation to infer the lens mass, 
imposing the effective surface threshold criterion for strong lens- 
ing to connect the lens mass to the halo mass. Finally, we used 
the milli-lensing optical depth in each scenario to calculate the 
expectation number of detected milli-lenses in the source sample 
of the SMILE project. 

We found that the probability of strong milli-lensing by DM 
haloes strongly depends on the model, being regulated by the 
properties of DM particles which dictate the inner structure of 
haloes, as well as their number density. We have shown that even 
within the CDM model, the lensing optical depth is quite sensi- 
tive to variations in the concentration-mass relation, leading to 
very different expectation values of detected milli-lenses in the 
SMILE source sample. Milli-lensing observations might there- 
fore enable us to constrain the concentration-mass relation down 
to sub-galactic mass scales. 

In addition, we have demonstrated that DM scenarios which 
are in favor of core-like density profiles, such as the SIDM 
one investigated here, are unlikely to produce milli-lenses be- 
cause they predict haloes with low-density centers. However, our 
method allows to also probe scenarios like core collapse which 
enhance considerably the probability of milli-lensing. 

Finally, we have shown that haloes consisting of WDM 
lead to an extremely small milli-lensing optical depth due to 
their combination of low concentration and mass-function cut- 
off. Even if a steeper concentration-mass relation (such as Eq. 
39) is used, the cutoff in the number density below sub-galactic 
scales still prevent the milli-lensing optical depth from increas- 
ing significantly. Therefore, the detection of milli-lenses would 
provide definitive evidence against the WDM model and more 
generally models that exhibit a cutoff in their halo mass function 
affecting the 10° — 10° Mo mass scales. 

In Fig. 6, we summarize our results, plotting the expectation 
value of detected milli-lenses for all models investigated in this 
study. The blue point corresponds to the CDM model A (Sect. 
4.1.1), while the green point refers to the CDM model B (Sect. 
4.1.2). The black point corresponds to the core collapse SIDM 
scenario C (Sect. 4.2.1) and the red point to the WDM model 
D (Sect. 4.3). Even among the limited number of DM models 
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Fig. 6. Expectation number of detected milli-lenses in the SMILE 
project. The error-bars have been derived assuming Poisson-like error 
in the calculation of the expectation values. 


studied here, milli-lensing observations of source samples com- 
parable to that of SMILE hold significant discriminating power. 

Even though sub-galactic haloes are expected to be almost 
empty of baryons, one source of uncertainty in our work might 
be the fact that we ignore the overall effect of baryons in the 
internal structure of haloes which in principle can alter the abil- 
ity of those haloes to act as milli-lenses. Although such systems 
are DM dominated, the baryons might play a crucial role in the 
strong lensing and therefore it is left for a future investigation, 
since the purpose of this paper is mostly to highlight the point 
that milli-lensing observations can be used to constrain the na- 
ture of dark matter and further discriminate between currently 
viable models. 

Another possible source of uncertainty arises from the fact 
that some DM haloes might host a SMBH at their center, which 
would contribute significantly in the strong lensing signal from 
those haloes, thereby modifying our results. However, neither 
the fraction of haloes that host such objects at their centers nor 
the accurate relation between SMBH and halo mass are known. 
Hence we do not take this possibility into account in this study, 
but we plan to explore this scenario in the future. 

Here, we assumed a certain redshift distribution for the ~ 1/3 
of sources for which no redshift measurements were currently 
available. This choice partially affects the results, although not 
our qualitative conclusions. For instance, if we instead assumed 
that missing redshifts are due to weakest and farthest sources, 
then the real redshift distribution for these sources would extend 
toward higher redshifts which correspond to higher values of the 
optical depth, thereby increasing the expectation number of de- 
tected milli-lenses and enhancing the probability of detecting a 
milli-lens. 

In this study, we have concluded that: 1) the source sample 
included in the SMILE project is sufficiently large to enable in- 
ferences about the nature of DM; 2) WDM haloes are highly un- 
likely to produce even a single strong milli-lensing event in the 
source sample of SMILE; 3) SIDM haloes can only act as strong 
milli-lenses in the case where self-interactions trigger the core 
collapse mechanism, leading to highly dense cores; and 4) the 
ability of CDM sub-galactic haloes to act as milli-lenses strongly 
depends on the mass-concentration relation. Finally, we have 
shown that if CDM is indeed the relevant model for describing 
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the properties of DM particles, then milli-lensing observations 
will enable us to further constrain the relationship between con- 
centration and halo mass down to sub-galactic mass scales. 
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Appendix A: Press-Schechter mass function 


In this appendix, we offer a short prescription for calculating the 
halo mass function. For more details on the subject, see Mo et al. 
(2010). 

Using the Press-Schechter formalism (Press & Schechter 
1974), the differential halo mass function is expressed as 


dou |-#2 


aN _ |2 Pm) ôk) a 
dM | P 2 


n(M,z) = — = 
( ) dM m M ow, 


| (A.1) 


where pPm(Z) = Qn(1 + z), 6-(z) is the overdensity of a struc- 
ture collapsing at redshift z linearly extrapolated to the present 
epoch, and oy accounts for the linear rms fluctuation (variance) 
of the density field on scale M. Under the spherical collapse as- 
sumption and concordance cosmology, the critical overdensity is 
given by (e.g., Eq. C30 in Pavlidou & Fields 2005) 


1 
6-(Z) = Do (A.2) 


where D(z) is the normalized linear growth factor, i.e., D(z = 
0) = 1, calculated by 


DQ) = G |20)" a +2)]/E[2)'?], (A.3) 
while w = Q4 /Qn and 
(2 Bs u>yt/2 u ip? 


The variance oy, normalized to be equal to og when R = 
847! Mpc, is given by (see Eq. 10 in Pavlidou & Fields 2005) 


ee fi dk POW? (KRMK 
2 = 


8 eS (A.5) 
J, dk P(k)W?(k8h-! Mpo)k? 


where the og parameter is a direct observable quantity, while 
W(kR) refers to the window function (filter) and R is the charac- 
teristic radius of the filter related to the mass M through 

1/3 

M 
ran =(=) 
YfPm 

where yy is a parameter depending on the shape of the filter. 
Here, we employ a sharp in k-space filter given by 


(A.6) 


pf 1 if|kl<R7 
Wek) = { Too E (A.7) 


Given this choice, yr becomes equal to 67? (Mo et al. 2010). 

In Eq. (A.5), P(k) stands for the linear matter power spec- 
trum. Invoking linear theory, we can write the matter power spec- 
trum to be 


P(k) œ T?(k)Pinit(k), (A.8) 


where P;,,;; is the initial power spectrum proportional to k” with 
n = 0.97, while T(k) corresponds to the transfer function. Bond 
& Efstathiou (1984) offer the following simple numerical for- 
mula for the transfer function T(k) which is consistent with the 
ACDM model (see also Jenkins et al. 2001) 


1 
[1 + [aq + (ba)??? + et] 


where q = k/T, T = Qyoh, v = 1.13, a = 64 ho! Mpc, b = 
3 h! Mpc, and c = 1.7 h~! Mpc. 


T(k) = (A.9) 
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